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Comparison of Hadronic Interaction Models at Auger Energies 

D. Heck a *, M. Risse a , and J. Knapp b 

a Institut fiir Kernphysik, Forschungszentrum Karlsruhe, P.O. Box 3640, D-76021 Karlsruhe, Germany 
b Dept. of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, United Kingdom 

The three hadronic interaction models DPMJET 2.55, QGSJET 01, and SIBYLL 2.1, implemented in the air 
shower simulation program CORSIKA, are compared in the energy range of interest for the Pierre Auger ex- 
periment. The model dependence of relevant quantities in individual hadronic interactions and air showers is 
investigated. 



1. INTRODUCTION 

The engineering array of the southern Pierre 
Auger Observatory (PAO) jj) in Argentina has 
recently started to take data. The interpretation 
of the measurements requires reliable numerical 
simulations of extensive air showers (EAS) in the 
extremely-high energy (EHE) range > 10 19 eV. 
A large uncertainty in such simulations arises 
from the models which describe the hadronic in- 
teractions. As one has to extrapolate in energy 
from accelerators by several orders of magnitude 
and into the forward kincmatical range which is 
unobserved by collider experiments, the models 
have to rely on theoretical guidelines to describe 
the EHE collisions. In this contribution model 
predictions for single interactions of p-p, p- 14 N, 
and 7r- 14 N collisions are compared, and we exam- 
ine how their features influence measurable EAS 
quantities such as the longitudinal development 
with the shower maximum A max and the lateral 
distributions of particle densities at ground. 
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Figure 1. Inelastic proton-air cross-sections of 
the considered models. Included are experimental 
data of air shower measurements. Further details 
see ref. Ol. 



by its authors. But already at this energy some 
interesting trends of this model show up. Fig. [j] 
shows the proton-air cross-sections for production 
of secondary particles as function of energy. 



2. MODELS 

Presently three hadronic interaction codes 
coupled with the EAS simulation program 
CORSIKA § are able to treat hadronic colli- 
sions at energies > 10 19 eV: DPMJET 2.55 [§Jj, 
QGSJET 01 [OTh and SIBYLL 2.1 §§. In this 
study neXus 2 [| is partly included despite of 
its upper limit of < 2 x 10 17 eV recommended 
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3. SINGLE INTERACTIONS 

All models are tuned to reproduce the avail- 
able collider data. They agree fairly well in the 
pseudo-rapidity, multiplicity, and transverse mo- 
mentum distributions obtained by the UA5 [Q, 
CDF [0, and P238 @ experiments. But when 
extrapolating to higher energies, already at £q a b 
= 10 17 eV (corresponding to y/s « 14 TeV, which 
will be reached by the future LHC-collider), the 
mid-rapidity density at \rj\ < 4 of QGSJET ex- 
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Figure 2. Pseudo-rapidity distribution of charged 
particles in p-p collisions at -E^b = 10 17 eV. 
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Figure 3. Distribution of charged particle multi- 
plicity in 7r- 14 N collisions at £i a b = 10 17 eV. 



ceeds that of neXus by up to a factor of «2 as 
shown in Fig. |^. In Fig. the charged particle 
multiplicity distribution of 7r- 14 N collisions shows 
a similar behaviour. Again QGSJET predictions 
are higher than those of all other models, while 
the neXus distribution repeats the low multiplic- 
ity already shown in the p-p collisions of Fig. ^. 
This behaviour reflects the different treatment 
of the Pomeron exchange by the various mod- 
els. In QGSJET a hard Pomeron is always cou- 
pled to the partons of projectile and target via 
soft Pomerons, which finally produce the large 
number of secondary particles. Secondaries pro- 
duced by cutting soft Pomerons appear in the 
mid-rapidity range and take away only a minute 
energy fraction. Therefore they influence the de- 
velopment of charged particle numbers in EAS 
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Figure 4. Distribution of longitudinal momen- 
tum fraction carried away by the most ener- 
getic baryon emerging from p- 14 N collisions at 
Slab = 10 19 eV. 

rather insignificantly. 

Already in the calculation of the Pomeron ex- 
change probabilities neXus applies strict energy 
conservation, which reduces the number of cut 
hard and soft Pomerons, thus reducing the over- 
all multiplicity. This decreases the number of par- 
ticles emitted in the pseudo-rapidity range below 
\r)\ < 6 while the fraction of secondaries produced 
in the very forward direction \r]\ > 8 resembles 
that of the other considered models. 
The distribution of the longitudinal momentum 
fraction taken away by the leading particle is dis- 
played in Fig. [TJ Large differences up to a factor 
of 10 show up in the diffractive region x\ a b > 0.8, 
but also around x\ a b ~ 0.4 significant deviations 
are visible. The extremely low probability of 
DPMJET to produce baryons at xuh « 0.97 in p- 
14 N collisions is remarkable. It reflects the insuf- 
ficient knowledge on high- and low-mass diffrac- 
tion: Any detailed treatment within the interac- 
tion codes is strongly model dependent. 

4. INFLUENCE ON EAS FEATURES 

In hadronic EAS the ir^ -mesons are the most fre- 
quent secondary hadronic particles and nitrogen 
is the most abundant component of air. Therefore 
the character of an EAS is essentially influenced 
by the features of 7r- 14 N collisions. In such inter- 
actions QGSJET exhibits a lower average elas- 
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Figure 5. Average of longitudinal momentum 
fraction carried away by the most energetic me- 
son emerging from 7r- 14 N collisions as function of 
energy. 

ticity than DPMJET and SIBYLL, see Fig. | 
Therefore EAS simulated with QGSJET develop 
faster than those produced with the other models. 
The second important influence on the develop- 
ment process arises from the cross-sections. The 
7r-air production cross-sections behave similar to 
those shown in Fig. 0, but with DPMJET stay- 
ing well below QGSJET above 10 15 eV. Therefore 
EAS simulated with the DPMJET model are ex- 
pected to develop late, while SIBYLL showers will 
fall somewhere in between QGSJET and DPM- 
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Figure 6. Longitudinal development of charged 
particle number. Vertical incidence, i?o=10 19 eV, 
E c > 0.1 MeV. The vertical dotted line indicates 
the vertical depth of the PAO in Argentina 0. 




Figure 7. Depth X max of shower maximum as 
function of primary energy. Vertical incidence, 
E =10 19 eV, E c > 0.1 MeV. References to exper- 
iments (fluorescence = full symbols, Cherenkov 
techniques = open symbols) are given in Ell . 

JET. This is demonstrated in Fig. || for proton 
and Iron induced showers of 10 19 eV. At the upper 
margin of Fig. ^ the positions of the maxima are 
indicated. They nearly coincide for Iron induced 
showers, but for proton induced EAS the differ- 
ence between QGSJET and DPMJET amounts 
to « 40 g/cm 2 . 

The depth of shower maximum AT max as function 
of primary energy is plotted in Fig. [?]. While 
for Iron induced showers even at highest energies 
the V max values coincide within < 20 g/cm 2 , 
with a tendency of DPMJET to predict deeper 
penetration, the lines of proton induced show- 
ers show a clear divergence with increasing en- 
ergy reaching 55 g/cm 2 at 10 20 eV. Especially 
the QGSJET distribution flattens with increas- 
ing energy. Starting at 10 15 eV with 65 g/cm 2 
per energy decade the slope reaches 45 g/cm 2 per 
decade at 10 20 eV. The flattening has to be at- 
tributed to the lowering of the 7r- 14 N elasticity 
(Fig. ||) together with an increasing 7r- 14 N cross- 
section of QGSJET which exceeds that of DPM- 
JET by 10% at E lah = 10 17 eV. Above 10 15 eV 
the V max slopes of the other models show no sig- 
nificant change for proton induced EAS. 
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Figure 8. Ratios of lateral particle densities rela- 
tive to QGSJET. Vertical incidence, £ =10 19 eV, 
ethin = 10~ 6 (optimized ||), E c > 0.25 MeV, 
> 0.1 GeV. 



Because of its early development QGSJET ex- 
hibits the flattest lateral distribution for electrons 
and for muons. This results in the highest particle 
densities at distances > 300 m not only for pro- 
ton but also for Iron induced EAS. Therefore we 
compare in Fig. |^ the lateral distributions of the 
other models relative to those of QGSJET. With 
increasing distance from the shower core the den- 
sities predicted by the other models are always 
smaller. The effect is more pronounced for proton 
induced showers. In general the SIBYLL densi- 
ties are lowest and the differences are larger for 
the muon distributions. 

5. CONCLUSIONS 

We emphasize the importance of future collider 
experiments to measure the very forward range 
of particle production thus reducing the present 
uncertainties in EAS simulations caused by the 
hadronic interaction models. Despite the un- 
known elemental composition of the primary cos- 
mic radiation at EHE, the PAO should be able to 
constrain the hadronic interaction models by eval- 
uating hybrid EAS events which combine the lon- 
gitudinal information measured by fluorescence 
telescopes and the lateral particle densities reg- 



istered by water Cherenkov surface detectors. 
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